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On the Unique Reconstruction of Induced Spherical 

Magnetizations 

Christian Gerhard^] 

August 4, 2016 

Abstract. Recovering spherical magnetizations m from magnetic field data in the exterior is a 
highly non-unique problem. A spherical Hardy-Hodge decomposition supplies information on what 
contributions of the magnetization m are recoverable but it does not supply geophysically suitable 
constraints on m that would guarantee uniqueness for the entire magnetization. In this paper, we 
focus on the case of induced spherical magnetizations and show that uniqueness is guaranteed if 
one assumes that the magnetization is compactly supported on the sphere. The results are based 
on ideas presented in [2] for the planar setting. 
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1 Introduction 

The lithospheric contribution to the Earth’s magnetic field is due to magnetized rocks in the 
Earth’s crust, which can be expressed as a spherical shell Q r ,R = {x G M 3 : r < |x| < R}. 
The magnetic potential V that is generated by a vectorial magnetization m : £l r ,R —> K 3 can 
be expressed by 


V(x) 


1 

47T 



x-y 
x - y I 3 


d y, 


x € M 3 . 


( 1 . 1 ) 


Recovering m from knowledge of V only in the exterior = {.x E M 3 : |x| > i?} is a 
highly non-unique problem. The relation (1.1) actually seems to be a vectorial version of the 
gravimetry problem (see, e.g., ummm) and reveals similar uniqueness issues. However, 
opposed to the gravimetry problem, where the assumption of a harmonic mass density leads 
to uniqueness, the assumption of a harmonic magnetization m would still maintain a certain 
non-uniqueness. The non-uniqueness even persists if we restrict ourselves to induced mag¬ 
netizations of the form m = Qv, where v denotes a known inducing vector field and Q an 
unknown scalar susceptibility. In |24| . it has been shown that a constant susceptibility in the 
spherical shell $l r ,R produces no magnetic effect in the exterior 0,^ * (for any inducing vector 
field of the form v = VC/, where U is harmonic in These considerations have been 

generalized to ellipsoidal shells in ns- A discussion of further examples of magnetizations 
and related uniqueness issues can be found, e.g., in [5j- 

Since the thickness of the spherical shell Q r _ r. where magnetization in the Earth’s litho¬ 
sphere occurs, is only a few tens of kilometers (thus, negligibly small compared to the Earth’s 
radius), it is geophysically reasonable to reduce the considerations to vertically integrated 
magnetizations m : Qr -a M 3 on the sphere f Ir = {x £ M 3 : |x| = R}. Therefore, from now 
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on, we consider the relation 


V(x) = 


4tt 


m(y) 


x-y 
\x - y\ c 


d u(y), x E 


( 1 . 2 ) 


instead of (1.1). By dw we denote the surface element on the sphere Qr. The non-uniqueness 
of recovering a vertically integrated magnetization m from the knowledge of V in can be 
characterized by a fairly well-known decomposition (see, e.g., dUHununi uni mi m i2i) 


m = fh^ 


fh^ , 


(1.3) 


which has the property that V = 0 in if and only if fh^ = 0 (in other words, any 
magnetization of the form m = fh W + fh® produces no magnetic potential in the exterior 
n&jrf). yy e ca jj suc ]- 1 a decomposition a Hardy-Hodge decomposition (cf. [3j for its Euclidean 
counterpart in M 2 ) and treat it in more detail later on. An illustration of this decomposition 


for recent magnetization models is supplied in [13] ■ Nonetheless, a characterization of m by 
(1.3) still states that the contributions fh ® and fh ® cannot be reconstructed from knowledge 


of V in fl'p’ 4 . Even if we assume an induced magnetization m = Qv, non-trivial susceptibilities 


Q have been constructed in m that generate a magnetic potential V via ( |1.2[ ) which vanishes 
in Qft* (again, assuming that v is known and of the form v = Vt/, where U is harmonic in 
underlining the non-uniqueness of the problem. 

In this paper, we show that an induced magnetization m = Qv is uniquely recoverable 


from (1.2) if V is known in fl??* and if one imposes the additional condition that m has 


compact support in Qr (it is not necessary that v is of the form v = Vt/). Now, if there 
exists a model of the vertically integrated induced magnetization that is very accurate in 
some region T/j C fl/j of the Earth, the residuum of the magnetic potential generated by 
this model magnetization and the magnetic potential obtained from actual global magnetic 
field measurements forms a magnetic potential that can be regarded as being generated by 
a magnetization with compact support in fl/j \ Tr. The residual magnetization can then be 
determined uniquely due to our result and, together with the accurate model magnetization 
in T/j, we obtain a trustworthy model for the induced vertically integrated magnetization in 
Qr. The proof of our uniqueness result is based on a combination of the spherical Hardy- 
Hodge decomposition from isi muni and the ideas presented in [2] for the case of thin-plate 
magnetizations in the plane M 2 . 

After introducing some notations in Section [l.l[ we briefly recapitulate results from [2tfl6j 
in Section 1.2 in order to highlight their relations to the spherical case later on. In Section[2j we 
introduce the spherical Hardy-Hodge decomposition and in Section [3j we discuss constraints 
on induced spherical magnetizations m = Qv that guarantee uniqueness if V is known, e.g., 
only in the exterior Furthermore, we discuss the numerical reconstruction of m and 

supply some examples in Section |4j 


1.1 Notations 

For the spherical setting, we assume that the magnetization m is located on the unit sphere 
Q = Qi and we denote the exterior by VL ext = {x E M 3 : |x| > 1} and the interior by 
= {x E M 3 : \x\ < 1}. In the Euclidean setting of thin-plate magnetizations in [Jj, the 
magnetization m is restricted to M 2 ~ M 2 x {0} and the “exterior” is represented by the upper 
half-space and the “interior” by the lower half-space M 3 . In the following, we essentially 
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use identical notations for the spherical setting and for the Euclidean thin-plate setting (of 
which the thin-plate setting only appears in Section 1.2). Vector fields mapping fl or M 2 
into M 3 are denoted by lower case letters f,g ,..., and if they are square-integrable, we say 
that they are of class l 2 (fl) or 1 2 (M 2 ), respectively. Scalar fields mapping Q or M 2 into M are 
denoted by upper case letters F,G,..., and if they are square-integrable, we say that they 
are of class L 2 (fl) or L 2 (M 2 ), respectively. By Latin letters x,y,... we mean vectors in M 3 or 
M 2 , by Greek letters £, 77 ,... unit vectors in 1L 

The surface gradient V* is defined as V* = 0), for x = (aq, X 2 ) 6 R 2 ~ K 2 x {0}, 

in case of the Euclidean setting. In the spherical setting, it is defined by the connection 


= rV x — to the gradient V a 


m 


where r = 


and ( = n 6 !1. The surface curl 


gradient L* reads as L* = (0,0,1) A V* = ( — g§^,0) in the Euclidean setting and as 

Lj) = £ A V| in the spherical setting, where A denotes the vector product. For convenience, 
we introduce the following notation for the (tangential) differential operators above: 


)( 2 ) = V* 


= L* 


(1.4) 


These operators are complemented by the operator which is given by 

o (1) = (0,0,1) id (1.5) 

in the Euclidean setting and by 

o^=£id (1.6) 

for the spherical setting (id denotes the identity operator), o^ always points in normal 
direction with respect to M 2 ~ M 2 x {0} or Q, respectively. Last, we need the Beltrami 
operator A* which is defined by A* = V* • V* = L* • L*. In the Euclidean setting this means 
that A* = (gf ^) 2 + (gf^) 2 > f° r x = ( x i, x 2) £ R 2 - t 2 x {0}, and in the spherical setting 
that the connection A| = r 2 A x — J^r 2 J^ to the Laplace operator A^, in M 3 holds true, where 
r = |x| and ^ = A G Q. 


1.2 The Thin-Plate Case 

For the thin-plate setting, we assume that the magnetic potential V is generated by a vectorial 
magnetization m of class 1 2 (M 2 ). Then we can write 

vw = r [ m (y)- \ X xeM 3 . (1.7) 

4tt 7 R 2 |x-(y,0)| d 

We use the following definition in order to characterize m with respect to its effect on V. 

Definition 1.1. Two magnetizations m, fh G 1 2 (M 2 ) are called equivalent from above if the 
corresponding magnetic potentials V and V (given by & are equal in the upper half-space, 
i.e., if V = V in M^.. They are called equivalent from below if V and V are equal in the 
lower half-space, i.e., ifV = V in M 3 . A magnetization m is called silent from above if it is 
equivalent from above to fh = 0 and silent from below if it is equivalent from below to m = 0. 
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A decomposition of m that reflects this characterization is the so-called Hardy-Hodge 
decomposition (for details on the thin-plate case and all results mentioned in this section, the 
reader is referred to 0)- For that purpose, we require the following vectorial operators: 


o« = (tfi,i? 2 ,id), (1.8) 

d (2) = (-«i,-i?2,id), (1.9) 

o (3) = (-^2^1,0) , (1.10) 

where Ri, i = 1, 2, are the Riesz transforms 

Ri[F](x) = lim f F(y)^ —^dy, x = (xi,x 2 ) 6 M 2 , (1.11) 

2vr J R 2 \ Bs ( x ) \x - y | d 


of a scalar function F of class L 2 (M 2 ). We can now formulate the following theorem. 

Theorem 1.2 (Hardy-Hodge Decomposition). Any function f = (F X1 , F X2 , F X3 ) E 1 2 (M 2 ) can 
be decomposed into 

f = / (1) + / (2) + / (3) = d (1) [A] + d (2) [F 2 ] + o^[F 3 ], (1.12) 


with scalar functions F\, _F 2 , F$ given by 

Fi = - (— Ri[F xl ] — R 2 [F X2 ] + F X3 ), ( 1 - 13 ) 

F 2 = - (Ri[F Xl ] + R 2 [F X2 ] + F x 3 ), ( 1 - 14 ) 

F 3 = R 2 [F X1 ] — R\[F X2 }. ( 1 - 15 ) 

Spherical counterparts to the above theorem are introduced in more detail in Section [2j 
A helpful formal notation for the operators o o^ 2 ', 0 ^ that emphasizes the connection 
between the spherical and the Euclidean thin-plate case is given in the next remark. 


Remark 1.3. It is well-known that the Riesz transform Ri can formally also be expressed 
as Ri = V*(A*)- 2 , i = 1,2, where V* denotes the i-th component of the surface gradient 
V* from Section 0 With this notation at hand we can now reformulate the operators 


( 1 - 8 )—( 1 . 10 ): 


o (1) = o (1) +o (2) ( A *)"2 j (1.16) 

o (1) = o (1) - o (2) ( A *)"5 } (1.17) 

o( 3) = o (3) (A*)"5 . (1.18) 


Supported by the Hardy-Hodge decomposition from Theorem 1.2 [0] have derived sev¬ 


eral characterizations and uniqueness results under the constraint of unidirectionality and/or 
locally compact support on the magnetization m. We list those who relate to results for the 
spherical setting later on in Sections [2] and [3j 


Theorem 1.4. Let m E 1 2 (M 2 ) and m^ 2 \ 
following assertions hold true: 


be given as in Theorem 


1.2 


Then the 
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(a) The magnetization m^ is equivalent from above to m while m ^ is equivalent from 
below to m. 

(b) The magnetization m is silent from above if and only if fh ( 2 ) = 0 while m is silent from 
below if and only if mS 1 ' = 0 . 

(c) If supp(m) C r, for a region Tel 2 with T 7 ^ M 2 , then m is silent from above if and 
only if it is silent from below. 

Corollary 1.5 (Unidirectional Magnetizations). Let m £ 1 2 (M 2 ) be a non-tangential unidi¬ 
rectional magnetization, i.e., m = Qv for a scalar function Q £ L 2 (M 2 ) and a fixed direction 
v = (vi,V 2 ,vs) £ M 3 , V 3 / 0. Furthermore, let T C M 2 be a region with T / M 2 and 
supp(m) C T. 

Then m is equivalent from above to no other unidirectional magnetization m = Qv with 
supp(m) C T. Analogously, m is equivalent from below to no other unidirectional magnetiza¬ 
tion m = Qv with supp(m) C T. 


2 Spherical Decompositions 

From now on, we are strictly working in the spherical setting, i.e., we are investigating the 
magnetic potential V that is generated by a magnetization m £ 1 2 (P): 

V(x) = — [ m(iq) ■ £—x £ M 3 . (2.1) 

4 ?r Jn - t?| d 


2.1 Vector Spherical Harmonic Representation 


A spherical version of the Hardy-Hodge decomposition from Theorem 1.2 has been used in 


geomagnetic applications for quite some time in form of a decomposition of vector fields with 
respect to vector spherical harmonics £{., £\, £\ (see, e.g., [U H31 HS1 ESI E2] ) • These 
vector spherical harmonics can be defined via a suitable connection to the inner harmonics 
H'lffi and the outer harmonics (i.e., the harmonic extensions of scalar orthonormalized 
spherical harmonics Y n ^ into P m * and Ll ext , respectively). More precisely, 


-.(i) 


x— 

xg n ext 


Vni (0 = * lim V x H™l(x) 

5 x — 7 


n £ No, k = 1 ,. 

.. 2 n + 1 , 

^ £ P, 

( 2 . 2 ) 

n £ N, k = 1 ,.. 

. 2 n + 1 , 

^ £ P, 

(2.3) 


x—t£ 

Igfi ™ 1 

with normalization constants jin' = (n + l)(2n + 2), fiff = n(2n + 1). The vector spherical 
harmonics y ^ ' k , n £ N, k = 1,... 2n + 1, are chosen such that, together with (2.2) and (2.3), 
they form a complete orthonormal function systems in 1 2 (P) with respect to the inner product 
(•, These properties imply that a square-integrable vector field / of the form / = VU, 

where U is harmonic in Ll ext , can be expressed by 


00 2n+l 


f = J2 m)£, 


,(1) 

,,k 


(2.4) 


n =0 k= 1 
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on the sphere fh Analogously, a square-integrable field / of the form / = VC/, where U is 
harmonic in f2 mt , can be expressed by 


oo 2n+l 


/ = E </> ^S) i2 ( n )^S ( 2 - 5 ) 

71=1 fc=l 

on the sphere fh More details on the involved (vector) spherical harmonics can be found, 
e.g., in EQj. Most general, one can state the following theorem. 

Theorem 2.1 (Spherical Hardy-Hodge Decomposition I). Any function f G l 2 (fl) can be 
decomposed into 

f = /<■>+ p>+p> 

oo 2n+l 

= Z Z d)i 2 (n)yS + Z Z Cfc)i 2 (n)Ci + &i a (n)P^l 


( 2 . 6 ) 


oo 2n+l 


oo 2n+l 


i(3) 

,,k' 


n =0 fc=l 


n=l fc=l 


n=l fc=l 


Remark 2.2. In order to investigate the consequences of Theorem 2.1 for the uniqueness of 
the magnetization m and the magnetic potential V in (2.1), we first observe that, for x G £l ext 
and 77 E f2, 


X — Tj 


1 


1-rv = lim V y - - r . 

\x — rj\ 6 v ^int \ x ~ y\ 


(2.7) 


Clearly, , 1, is harmonic in and (2.5) implies a representation of the form 


x — Tj 
\x — ?/| s 


oo 2n+l 

EE 

n=l fc=l 


\X 


> Cfe) Ci fa) > * G next > v e 

'I / l 2 (fi) 


R 


( 2 . 8 ) 


In detail, using the addition theorem for scalar spherical harmonics and a series representation 
of 1 in terms of Legendre polynomials P n , we obtain 

\X -\ 


X — 71 „ 1 \ - 

- = Inn 

\o y—pn y \%\ 


\x — 7/1" 


yen 1 ' 


72—0 


\y\ 


Pn 


y_ 

\y\ 


(2.9) 


47T 


2ji+1 


“55 E pp E w (9) 


2/^77 / 2n + 

y eQ. %rit 72—1 
OO 272+1 


fc=l 


EE 


/ 47r(/2 : 


(2)nI 


—■' ' \ 2 n + 1 

=i k= l \ 


h: 


ext 


i (=7)) v% w 


If we now substitute (2.8) or (2.9) into the representation (2.1) of V, we see that, due to the 
orthonormality of the vector spherical harmonics, the contribution rh^ of m = rh W + rh^ + 
?7T.( 3 ) generates the exact same magnetic potential V in Q ext as the entire magnetization m 
itself. Computations analogous to (2.7), (2.8), and (2.9) for the case x G and jjGSI lead 
to the conclusion that the contribution inP of m = rh W + ih^ + generates the exact 
same magnetic potential V in f V nt as the entire magnetization m itself. 


From Remark 2.2 it becomes clear that the spherical decomposition of T heorern 1 2.1 1 re veals 
the desired properties corresponding to the thin-plate case in Theorem |1.4[ a),(b) (which has 
already been observed in P3). 
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2.2 Operator Representation 


In order to be able to o btain a spherical version of Theorem |1.4[ c), we reformulate the 
decomposition of Theorem 2.1 in terms of a set of pseudo-differential operators o^ 1 ), d^ 2 \ 5^ 


as indicated in mnDununi. More precisely, 


odW 1 ) (d+II-oW, 

5< 2 >=o< 1 ) (d-T+oW, 

a< 3 > = o (3 t ' 

where D denotes the pseudo-differential operator 



( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


(2.13) 


The previously introduced vector spherical harmonics can then be alternatively expressed by 

Vnl = n G No, k = 1,... 2n + 1, (2.14) 

Vni = (fin ) ) _ ^o (2) F n)jt , n <E N, k = 1,.. ,2n + 1, (2.15) 

Vnl = (^ 3} )-^ (3 )y n>fe , n € N, k = 1,.. .2n + 1, (2.16) 

with normalization constants pi}" 1 = (n + l)(2n + 2), jjffl = n(2n + 1), and pffl = n{n + 1). 
In other words, the properties of the decomposition from Theorem |2.1| carry over to a de¬ 
composition with respect to the operators oW, d^ 2 \ 5^. To formulate such a decomposition, 
we first recapitulate the spherical Helmholtz decomposition (see, e.g., Elina and references 
therein). 

Theorem 2.3 (Spherical Helmholtz Decomposition). Any function f 6 1 2 (D) can be decom¬ 
posed into 

f = + o (2) [F 2 ] + o (3) [F 3 ], (2.17) 


where the scalar functions F\, F 2 , F 3 are uniquely determined by the conditions F 2 (r/)dca(r/) = 

fn F 3(ri)du{rj) = 0 . 

The spherical Helmholtz decomposition simply states a decomposition of a spherical vec¬ 
tor field into a normal component and two tangential components (of which one is surface 
divergence-free and the other one surface curl-free). From |81 llllfT2] , we now take the following 
theorem. 


Theorem 2.4 (Spherical Hardy-Hodge Decomposition II). Any function f £ 1 2 (D) can be 
decomposed into 


f = fi + h + h = d (1) [Ti] + oW[F 2 ] + o^[F 3 ], (2.18) 
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where the scalar functions F\, F^, F 3 are uniquely determined by the conditions f n Fi(rj) — 
F 2 (rj)duj(r /) = J n F 3 (r])dLo(r]) = 0. If F\, F 2 , F 3 are the Helmholtz scalars of f as given in 


in 

Theorem\2.S\ then 


F 1 = -(B-^[F 1 ]-F2 + -B- 1 [F, 


F 2 = - ^D- i [ J F 1 ] + F 2 + -D“ i [ J F, 
Hi = F 3 . 


(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 


A slight modification of the operators b^ 1 ), 5 ^, o® highlights the relation of Theorem 


2.4 


to the Euclidean thin-plate case in Theorem |L2| and Remark |1.3[ Application of the operator 
(D +|) to and application of (D — |) to 0 ^ and leads to the operators 


o« = - 0 ( 2 ) 


o( 2 ) = 


( 2 ) 


= 0^> + O 


o( 3 ) = 0 ( 3 ) ( D -I 


D+- 

2 

D —- 
2 

-1 


-1 


( 2 . 22 ) 

(2.23) 

(2.24) 


-1 


Comparing this to (1.8)—(1.10) and (1.16)—(1.18), we see that o^ (D+j) and o^ (D—f) 
respectively, take the role of the Riesz transform in the thin-plate case. However, it should 
be remarked that in the literature the Riesz transform on the sphere is typically given by 
c>( 2 )(A *) - 2 (see, e.g., [6]). Furthermore, (D — is well-defined only if restricted to the 
space Lq(H) = {F E L 2 (H) : F(rj)duj(r]) = 0} since the constant functions form the 

nullspace of D — \ Yet, this is no restriction for our further considerations. As an alternative 
to Theorem 1 2.4 [ we can now state the following theorem. 


Theorem 2.5 (Spherical Hardy-Hodge Decomposition III). Any function f E l 2 (fl) can be 
decomposed into 

f = h + h + h = o (1) [A] + d (2) [F 2 ] + d (3 )[F 3 ], (2.25) 

where the scalar functions F\, F 2 , F 3 are uniquely determined by the conditions f n F 2 (rj)duj(r]) = 
F 3 (rj)duj(?j) = 0. If F\, F 2 , F 3 are the Helmholtz scalars of f as given in Theorem 
then 


2.3, 


1 


1 


1 


F 1 = -\F 1 + - D _ 1 [Fi] - D[F 2 ] + - 


F 2 = x Fi - - D _i [Fi] + D[F 2 ] - - D" 1 ^ 


f 3 = D[F 3 ] - 2 ^ 3 . 


(2.26) 

(2.27) 

(2.28) 


Proof. For F\ = (D+ 2 )[Fi], ^2 = (D — \)[F?\, and F 3 = (D — 2 )[F 3 ], we directly obtain 
the representations (2.26)—(2.28) from the corresponding representations of F\ , F 2 , and F 3 
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in Theorem 2.4 Concerning the uniqueness, we can restrict the considerations to the case 


/ = 0. Using fl2.22D - fl2.23D in fl2.25P leads to 
0 = d ( 1 ) [Fi] + d ( 2 ) [F 2 ] + b ( 3 ) [F 3 ] 


= 0 « (Fi + F 2 ) + oW (d +0 [A] + (° -0 [U]j + o (3) T 3 - 


(2.29) 

(2.30) 


The uniqueness of the normal part of / leads to F\ = — F 2 . Thus, the assumption f n F 2 (rj)duj(r /) 
0 implies Fi(r/)dui(r/) = 0 and, consequently, 


-o-l D + 2 


-1 


[A](r?)+ D- 


[F 2 ](rj) du(r]) = 0. 


(2.31) 


The uniqueness of the spherical Helmholtz decomposition in Theorem 2.3 then leads to 
(D+l) _ 1 [Fi] + (D-l) _ 1 [Fi] = (D+1) _ 1 [A] - (D-1) _ 1 [A] = 0 , or in other words, by 
application of D — 1, 


2-ID+i 


-1 N 


[Fi] = 0. 


(2.32) 


Since 2 — (D+^) 1 is injective on Lq(O), we get F\ = —F 2 = 0. Moreover, the condition 
fn U^dcu^?) = 0 and the uniqueness of the Helmholtz scalars in Theorem |2.3| implies F3 = 0 
and concludes the proof. □ 


Remark 2.6. We need to emphasize that the functions / 1 , f 2 , / 3 in Theorems 2.1 and 2.4 


and the functions f\, f 2l / 3 in Theorem 2.5 are identical. The theorems only differ in the 
applied operators and the corresponding scalar functions F \, F 2 , T 3 , and F \, F 2 , F 3 . For our 
later considerations, we work with the operators d^ 2 \ o® and the functions F\, F 2 , F 3 . 


3 Uniqueness Issues for Spherical Magnetizations 

We start by defining the notion of equivalence of magnetizations, analogous to the thin-plate 
case in Definition o 

Definition 3.1. Two magnetizations m, fh £ 1 2 (0) are called equivalent from inside if the 
corresponding magnetic potentials V and V (given by ( 2 . 1 )) are equal in the interior Ll int , 


i.e., if V = V in Ll int . They are called equivalent from outside if V and V are equal in the 
exterior Ll ext , i.e., if V = V in Ll ext . A magnetization m is called silent from inside if it is 
equivalent from inside to fh = 0 and silent from outside if it is equivalent from outside to 
fh = 0. 

Theorem |2.4| eventually allows the following characterizations of spherical magnetizations. 

Theorem 3.2. Let m E 1 2 (P) and m^ l \ m^ 2 \ fh^ be given as in Theorem 2.4 Then the 
following assertions hold true: 

(a) The magnetization mf 2 ^ is equivalent from outside to m while m^ 1 ’ is equivalent from 
inside to m. 
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(b) The magnetization m is silent from outside if and only if m^ 2 > = 0 while m is silent 
from inside if and only if = 0 . 

(c) If supp(m) C r, for a region Tc with T ^ O, then m is silent from outside if and 
only if it is silent from inside. 


Proof. Parts (a) and (b) are direct consequences of the considerations in Remark 2.2 Con¬ 
cerning part (c), we assume that m is silent from outside and can now conclude that nv 2 ^ = 0 
and, therefore, 


M 2 = i (V 1 ^] +M 2 + l - D ~ l [M 2 ] ) = 0, 


(3.1) 


where Mi, M 2 , M 3 are the Helmholtz scalars of m as supplied in Theorem 2.3 Plugging (3.1) 
into (2.19) implies M\ = — M 2 . In other words, 

m = fh*' 1 ) + = —o^M 2 + M 3 . (3.2) 

Expanding the Helmholtz scalar M 2 in terms of spherical harmonics leads to the expression 
o( 1 ) M 2 = Y^i{ M 2 ,yn,k)L 2 {Q.) 0 ^y n ,k- Next, for x G Q ext , we set 


00 2n+l 


U(x) = E E (M 2 ,Y n , k ) L 2 m H e n % (x) 

n =0 k =1 

p / 00 2n+l \ 

E E H nl (*) Y n,k (:n) M 2 (r,)Mv) 


(3.3) 


\n=0 k =1 
f 00 


in 

1 

47T 


E (n+i)p « • vj j M 2 (ri)duj(ri) 


\n=0 


1 — \x\ 


:M 2 (j])du>(ji). 


'n (1 4 - \x\ 2 — 2x ■ ?y ) 2 

The last equalities can be found, e.g., in m and references therein. Clearly, U is harmonic 
in H ext . Remembering (2.2), (2.3), and (2.14)—(2.16), we can conclude that 

(gH. 


lim V x U(x) = -4 1 } M 2 (0, 

X —>£ s 


(3.4) 


From (|3.1[) we find M 2 = — (D+|) 1 Mi, and because (D+^) 1 is selfadjoint, it holds 


i\-i 




2 > 1 ’ “ * 1 2^ 
-1 1 1_1 2 




1 — x 


(1 + \x\ 2 — 2 x ■ rj) 2 


Mi(r})duj(rj). 


(3.5) 


Since supp(m) CT, we get M\ = 0 in T c = H \ T and can extend U across T c to a harmonic 
function in the cone Cr c = {1 £ l 3 \ {0} : dh G r c }. 

Next, we observe that there must exist some function W on r c such that rh ^ = o^W 
in r c . This follows from (3.2) because m = 0 in T c , because rh^ is tangential, and because 
rh W is solely composed by a normal part due to o ^ and a tangential part due to o^ 2 \ Now, 
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setting W(x) = W{x/\x\), for x E Cr=, we see that W is harmonic in C r=. Combining our 
findings up to now, we have a function U -\-W that is harmonic in Crc and which satisfies 


m(£) = lim V x {U(x) + IU(x)) = 0, £ E T c . 

x£Q ex t 


(3.6) 


Therefore, the normal and the tangential derivative of U + W vanish on T c . Consequently, 
since U + W is harmonic in Cr c , typical analytical continuation arguments lead to U + W = 0 
in Cr=. In particular, ^ • \7 x {U(x) + W(x)) = A • V x U(x) = 0, for x E Cr=. Since VC/ 
is harmonic in 54 e:rf by the construction in (3.3), we obtain ■ \7 x U(x) = 0, for x E Ll ext . 
Combining this with (3.4) leads to 


lim ■ V x U{x) = -£ • o£ ) M 2 (Z) = -(d + 1 ) [M 2 ](0 = 0, 
*->•£ \x\ s V 2 / 

Egfiemt 11 


e e n. 


(3.7) 


In other words M 2 = 0. Observing (3.1), we now find D = 0, which leads to M\ = 0 

and eventually to = 0 by (2.19). Thus, 


m = fh^\ 


(3.8) 


which, by parts (a) and (b), implies that m is silent from outside as well as inside. 


□ 


Theorem 3.2 has some direct consequences for the uniqueness of induced magnetizations 
of the form 


m{Z) = Q(£, MO; ^ e 54, 


(3.9) 


where v is a known inducing vector field and Q an unknown scalar susceptibility. If supp(m)c 
r and if there exists another magnetization m = Qv with supp(m)c T that generates the same 
magnetic potential in Q ext , then m — m is silent from outside and Theorem |3.2| implies that 
m — m = o^M for some scalar function M. In particular, m — m = (Q — Q)v must be 
tangential. If the known inducing field v is non-tangential (i.e., u(£) • £ 7 / 0 for some £ E 54), 
this implies m = m. Put in more rigorous terms, we obtain the following corollary. 

Corollary 3.3 (Uniqueness of Induced Magnetizations). Let m E 1 2 (5I) be of the form m = 
Qv (where Q E L 2 (54) is unknown and v E l°°(fi) is known) with supp(m) C T for a fixed 
region T C 54 with T 7 / 52. Furthermore, we assume that v is non-tangential. 

Then m is equivalent from outside to no other induced magnetization m = Qv E l 2 (54) with 
supp(m) C T. Analogously, m is equivalent from inside to no other induced magnetization 
m = Qv E l 2 (54) with supp(m) C T. 


The corollary above suits an application to the induced Earth’s crustal magnetization 
since the inducing field v is typically the Earth’s core magnetic field, which is non-tangential 
and known from models such as [23j. A restriction to unidirectional magnetizations as in 


Corollary 1.5 (which demands further structure on the inducing field v but, in exchange, 
does not require to know v in advance) for the Euclidean thin-plate case has no particular 
relevance in a global spherical context. Actually, depending on how one defines spherical 
unidrectionality, a spherical counterpart to Corollary |1.5| does not necessarily hold true. 
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Figure 1: Example of vector fields m according to ( |3. 10 ), with Q = 1, ( = (0,0,1), and 
v\ = V 2 = V 3 = l (left) and with Q = 1, £ = (0,0,1), and v\ = V 2 = 0, V 3 = 1 (right). 


Remark 3.4. A fairly natural definition is to call a spherical magnetization unidirectional if 
it is constant with respect to a certain local spherical coordinate system. For simplicity, we 
regard magnetizations of the form 

m(0 = Q(0 (viC + V2(C -(C- OO + vsZ AC), (GO, (3.10) 

where C E 0 and ui, u 2 , V 3 E M are fixed and Q E L 2 (0) (an illustration is given in Figure [I]). 
Note that the tangential components £ — (C'OC an d £A£ vanish as £ approaches the poles ±£. 
Therefore, the vectors £, £ ~ (C • £)£, and £ A £ do not form a local coordinate system in the 
points £ = ±£. Yet, the example derived below can serve as a counterexample to uniqueness 
for spherical unidirectional magnetizations if we assume that the function Q vanishes around 

±C- 

A consequence of a spherical counterpart to Corollary 1 1.5 1 would be that the trivial mag¬ 
netization m = 0 is the only spherical unidirectional magnetization that is silent from outside 
and has support supp(m) C T. This particularly includes purely tangential magnetizations 

m(£) = Q(£)(v2(C-(C-00 + ^AC), (3.n) 

If such a magnetization is silent from outside and supp(m) C T, then the proof of Theorem 
Applying the Beltrami operator and observing the properties of M 2 leads to 


|3.2| implies that M 2 = 0, where M 2 is the Helmholtz scalar of m as indicated in Theorem 2.3 


0 = A*M 2 (£) = V* • m(£) (3.12) 

= v*q(£) • (v 2 (c - (£ • cm+vst a 0 + Qm {-2 v 2 (c • o+• (£ a 0) 

= C • (^V*Q(£) - V 3 V Q(C )) - 2u 2 (C • £)Q(£). 


We now choose £ = (0,0,1). Parametrizing the sphere 12 with respect to polar distance 
t E [—1,1] and longitude (p E [0, 2tt), and using representations of V* and L* with respect to 
this parametrization (see, e.g., HUE]), equation (|3.12[) can be rewritten as follows: 


d 


d 


0 = u 2 (l -t 2 )—Q(t,<p) - v 3 —Q(t,ip) - 2 v 2 Q(t,ip). 

We further simplify the above by choosing u 2 = 0 and obtain 

d 

0 = V3-^Q(t,ip), 


(3.13) 


(3.14) 
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implying Q(t,ip) = Q (t), for some function Q that only depends on the polar distance t. 
Summarizing the above considerations, we can conclude that a function Q E C^([— 1,1]) 
with Q (t) = 0, for t = £ • ( and £ £ T c = 11 \ T, yields a non-trivial spherical unidirectional 
magnetization 

m(0 = v 3 Q{£, ■ ()£, A C, £ G fi, (3.15) 

that is silent from outside and satisfies supp(m) C T. In other words, we cannot expect to 
obtain a spherical counterpart to Corollary |1.5| for spherical unidirectionality as meant above. 

On the other hand, if we call a spherical magnetization unidirectional if m(£) = Q{fi) v , for 
a fixed v = (vi,V 2 ,vs) E M 3 \ {0} (i.e., if we use an identical definition of unidirectionality as 
in the planar case of Corollary |1.5[ ), we were not able to confirm nor to disprove uniqueness. 

Remark 3.5. The existence of an induced magnetization m = Qv with supp(m) C T that 
is equivalent from outside to a (not necessarily induced) magnetization rn E 1 2 (H) with 
supp(m) C T is generally not guaranteed. Assume, e.g., that v(£) = (, ( E 1 1. If m and 
m are equivalent from outside, then m — m is silent from outside and supp(m — m) C T. The 
proof of Theorem |3.2| yields that this holds true if and only if M\ = Mi and M -2 = M 2 , where 
Mi, M 2 , and Mi, M 2 are the Helmholtz scalars of m and m, respectively. However, for our 
choice of v it is clear that M 2 = 0. In other words, for any magnetization m with M 2 ^ 0 and 
supp(m) C T there does not exist an induced magnetization m = Qv with supp(m) C T that 
is equivalent from outside to m. 

The situation of existence of an induced magnetization changes if we drop the condition 
that m = Qv has to satisfy supp(m) C T, at least if the inducing field v satisfies certain 
conditions. 


Definition 3.6. We call a vector field v E 1°°(H) admissible if |£ -u(£)| > C, £ E H, for some 
constant C > 0, and if the coefficients 

v n ,k,m,l= [ Y m ,i(r )) ^ • V*Y nk (rj)du(rj) (3.16) 

Jn V ■ v{ri) 

satisfy the following property: For every n 3 -summable sequence g n .k, nENo, fc = l,..., 2n+l, 
*- e -> Yl’^ > =oYl'k=i 1 n6 \9n,k\ 2 < OO; the infinite dimensional system of linear equations 



1 


1 

n+\ 



Tm,/ 9n,ki 


n E No, k = 1, ..., 2n + 1, 


(3.17) 


has a n 2 -summable solution 7 n ^, n E No, k = 1,..., 2n + 1, i.e., Ylk^i 1 n4 \ln,k\ 2 < 00 . 

Here, 5 U)m denotes the Kronecker delta. 

Corollary 3.7 (Existence of Induced Magnetizations). Let m E h 2 (fl) = {/ E 1 2 (H) : 
Si=iI]^oSfc=t ln 4 K/’2/i*fc)i 2 (^)| 2 < °°} be a given magnetization. Then, for every ad¬ 
missible vector field v E 1°°(H), there exists a Q E L 2 (H) such that the induced magnetiza¬ 
tion m = Qv is equivalent from outside to m. Analogously, for every admissible vector field 
v E 1°°(H), there exists a Q E L 2 (H) such that the induced magnetization m = Qv is equivalent 
from inside to m. 
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Proof. According to Theorems 3.2 and 2.5, m and m are equivalent from outside if and only 
if 


1 


D [Mi] + M 2 + - D [M 2 ] = 2 M 2 , 


(3.18) 


where Mi, M 2 are the Helmholtz scalars of m and M 2 the Hardy-Hodge scalar of m. Further¬ 
more, since v is admissible, we get Q(£) = j , £ £ Q, and, therefore, 


1 


1 


<M 2 ,y^) L » (n) = A*r„, t ) L2(n) = (3.19) 


/ 1 ^n,fc)l 2 (f2) 

n(n + 1) v ’ 

1 




n(n + 1 ) Ja V v(v) 

00 2m+l 


\ ^ ^ r 

I 1 \ E / ^m,l(v) 


n(n + 1) 


m =0 1=1 
00 2m+l 


r] ■ v{r /) 


V*Y n , k {rj)du(r,) 


n (n + l) ^ X/ 
' ' m =0 2=1 


Observing that the operator D 1 acts via 

00 2n+l 


D _1 [M 1 ] = E E —T( M i,^, fe ) L2( n)^ )fe 


n 1 1 U + 9 

n=0 fc=l 2 


(3.20) 


and using (3.19) in (3.18) leads to the following infinite dimensional system of linear equations 
for the Fourier coefficients (Mi, Y n ,k)i, 2 (n) of Mi: 

/* ' /* ' ( 1 . v n,k,m,l H ——[fin,m3k,l ) (Ml, ^m,j)L 2 (f2) = 2(1W2, ^n,fc)L 2 (f2)> (3.21) 

Ti\Ti o) ^ ■ 9 / 

n £ No, A: = 1,..., 2n + 1. 


m =0 2=1 


The admissibility conditions on v guarantee that this problem is solvable and that the ob¬ 
tained £)(£) = , £ £ O, lies in L 2 (D). The Fourier coefficients (M 2 , Y n ,k)i, 2 (n) °f M 2 


can be obtained via (3.19). Here, m £ h 2 (D) guarantees a sufficient decay of the Fourier 


coefficients (Mi, Y n ,k)i, 2 (n) suc h that the interchange of the series and integration in (3.19) is 
allowed. The resulting induced magnetization m = Qv is then equivalent from outside to m. □ 


Remark 3.8. If we choose u(£) = £, £ £ D (which served as a counterexample in Remark 


3.5 for the case supp(m) C T), the conditions on v of Definition 3.6 are clearly satisfied 


since v nik ,m,l = 0. In other words, for any magnetization m £ 1 2 (Q) there exists an induced 
magnetization m = Qv that is equivalent from outside. For general vector fields v it is more 


difficult to check the conditions from Definition 3.6 Here, the reader is, e.g., referred to [25 
and references therein. 
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Figure 2: The true susceptibility Q (left) and the inducing vector field v (right). 


alpha=0 



alpha = 10" s 



alpha = 10" 3 





Figure 3: The reconstructed susceptibilities Q for different values a = 0,10 6 ,10 3 (the 
images represent a frontal view of the sphere on which Q is plotted). 



4 Numerical Illustrations 


In this section, we want to illustrate that the theoretical result from Corollary 3.3 has an 
actual influence on the numerical reconstruction of induced magnetizations m = Qv. We 
assume to know the magnetic potential V on Qr (we choose R = 1.1 in this example) that is 
generated by the induced magnetization m = Qv, where 


Q(0 = 


(£-C) 4 , £-C<o, 

0, else, 


«(0 = £ + C - (C ■ OC. 


(4.1) 

(4.2) 


and ( = (0,0,1) is fixed, i.e., the magnetization m has compact support in the lower hemi¬ 
sphere r = {£ E : £ • £ < 0} (see Figure [2] for an illustration). 

In order to approximate the true (but unknown) magnetization m = Qv by some m = Qv, 
we denote by V[Q\ the magnetic potential that is generated by m and minimize the functional 


R[Q\ = \\ V [Q] - V\\ 1 2{ ^ r) + a||Q V ||f 2(nNr) . 


(4.3) 


The first term in (4.3) simply represents a data misfit that measures the deviation of V\Q\ from 


the known magnetic potential V. The second term in (4.3) penalizes magnetizations m = Qv 


that have contributions outside T, i.e., magnetizations that do not satisfy supp(m) C T. For 
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the numerical minimization of J-[Q\, we expand Q in terms of Abel-Poisson kernels (cf. [9j): 


N 


Q(0 = J2^ K (C^n) 

l-h 2 


n= 1 


K(H • Zn) = 


(1 + h 2 - 2h(£ ■ tn))2 


(4.4) 

(4.5) 


where h E (0,1) is a fixed parameter (influencing the localization of K ) and E 14, n = 
1,,N, are some predefined points indicating different centers of the kernel K (in our case, 
we choose N = 10235 and points that are uniformly distributed on 11). Under these 
conditions, the minimization of T\Q\ reduces to solving the set of linear equations 

M 7 = g, (4.6) 

where 


M — Q^ n ,rn)n,m=l,...,N 

= ([ V n (x)V rn (x)duj(x) + a [ \v(r])\ 2 K(ri ■ £ n )K(r] ■ £ m )duj(ri) 


T = ('Ym)ra=l,...,7V5 


/n \r 


n,ra=l,...,7V 


9 = {9n)n=l,...,N = ( / V n {x)V (x)du(x) 

\JQjr / n=l,...,7V 


(4.7) 

(4.8) 

(4.9) 


and 


Vn(x) = ■]- [ K(rj ■ U)v(v) ■ A — 73 doj(r)). (4.10) 

4 tt Jn \ x — V\ 6 

The quadrature rules from HI E| are used for the numerical evaluation of the occurring 
integrals. The reconstructed susceptibilities Q for different choices of a are shown in Figure 
[ 3 ] (the parameter h of the Abel-Poisson kernels is set to h = 0.9). One can see that for 
parameters a > 0, we obtain a fairly good reconstruction of the underlying true susceptibility 
Q, while for a = 0 (i.e., no penalization is taken into account for magnetizations m that 
violate supp(m) C T), we obtain an entirely different susceptibility Q. Latter generates the 
same magnetic potential on O/j as Q but it does not satisfy supp(rh) C T. 

In a second test example, we choose a slightly more complicated magnetization m = Qv. 


Q(0 


f 1000 (2 + £ • e) 4 cos(27r£ • C) sin(27r£ • C), £ • C < and £ • ( > 
\ 0, else, 

(C-Of+ C-(C-0£, 


(4.11) 

(4.12) 


where ( = (0,0,1) and ( = (0, \) are fixed. The magnetization m is again supported 

in the lower hemisphere r = {£Efl:£-£<0}, but the actual support is only a subset of T 
(see Figure [4] for an illustration). V denotes the magnetic potential on fl/j produced by m. In 
order to approximate the (unknown) magnetization m = Qv by a magnetization rh = Qv from 
knowledge of V, we again minimize J~[Q} as in (4.3)—(4.10). The parameter h of the involved 
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Figure 4: The true susceptibility Q (left) and the inducing vector field v (right). 



Figure 5: The reconstructed susceptibilities Q for different values a = 0,10 6 ,10 3 (the 
images represent a frontal view of the sphere on which Q is plotted). 


Abel-Poisson kernels is chosen to be h = 0.95 (in order to supply better localized kernels 
I\). The results are shown in Figure [5j Once more, for a = 0, we obtain a magnetization 
m that is equivalent to m from outside but does not satisfy supp(m) C T. Choosing a > 0 
leads to a reconstruction of the desired magnetization. However, we see (in particular in the 
center image of Figure [5]) that there are some undesired artefacts around the South Pole. 
Considering that these artefacts are reduced for an increased a, they might be due to the 
fact that the operator that maps the magnetization m on Q to its magnetic potential V on 
Qr has an unbounded inverse (the classical ill-posedness that typically occurs for potential 
field problems). In this paper, our focus is on ill-posedness in the sense of non-uniqueness. 
Some more sophisticated regularization methods that also address the ill-posedness due to 
unboundedness of the inverse for the planar magnetization problem can be found, e.g., in [Si- 
Last, it should be noted that in the second example, the lower hemisphere T that we chose 
for the numerical reconstruction is significantly larger than the actual support of m. This 
shows that the numerical scheme also works if we do not know the exact support of m. 

5 Conclusion 

We proved that for induced spherical magnetizations (where the inducing vector field is 
known) the additional assumption of compact support in some region T C H yields uniqueness 
for m. The numerical examples indicate that including this additional condition in the re¬ 
construction procedure guarantees picking the ’correct’ magnetization out of those that could 
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generate the measured magnetic potential. 
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